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Abstract: The full set of cosmological observables coming from linear scalar and 
tensor perturbations of loop quantum cosmology is computed in the presence of 
inverse-volume corrections. Background inflationary solutions are found at linear 
order in the quantum corrections; depending on the values of quantization parame- 
ters, they obey an exact or perturbed power-law expansion in conformal time. The 
comoving curvature perturbation is shown to be conserved at large scales, just as in 
the classical case. Its associated Mukhanov equation is obtained and solved. Com- 
bined with the results for tensor modes, this yields the scalar and tensor indices, 
their running, and the tensor-to-scalar ratio, which are all first order in the quantum 
correction. The latter could be sizable in phenomenological scenarios. Contrary to 
a pure minisuperspace parametrization, the lattice refinement parametrization is in 
agreement with both anomaly cancellation and our results on background solutions 
and linear perturbations. The issue of the choice of parametrization is also discussed 
in relation with a possible superluminal propagation of perturbative modes, and 
conclusions for quantum spacetime structure are drawn. 
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1. Introduction 

Loop quantum cosmology (LQC) [|J] provides the framework for implementing sev- 
eral effects seen to arise for the quantum geometry of loop quantum gravity (LQG) 
[0, 0, fl in a cosmological setting. Strict formulations of loop quantum cosmology are 
defined on minisuperspace, where one quantizes homogeneous spacetimes using the 
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methods of loop quantum gravity. The characteristic effects of discrete spatial geom- 
etry then emerge also in the reduced context, changing the dynamics of expanding 
universe models. The dynamics changes in particular at high densities, giving rise 
to mechanisms avoiding classical singularities. 

Quantum corrections to the classical spacetime structure not only exist at high 
densities but remain present, in weaker form, as the universe expands and dilutes. In 
such a regime, non-linearities as well as corrections from loop quantum gravity can 
be treated perturbatively in a gauge-invariant way. This has been done for linear per- 
turbations around spatially flat Friedmann-Robert son- Walker (FRW) models with 
one particular class of corrections expected from loop quantum gravity. These cor- 
rections, related to spatial discreteness, arise whenever an inverse of a certain metric 
component (or rather, densitized-triad components) appears in a matter Hamilto- 
nian or in the Hamiltonian constraint of gravity. Such inverse components are ubiq- 
uitous, for instance in kinetic matter terms, and thus the resulting corrections, called 
inverse- volume corrections, are an unavoidable consequence of loop quantum gravity. 
Specifically, the corrections are due to the fact that the quantized densitized triad 
has a discrete spectrum, with the value zero contained in the spectrum. Such an 
operator does not allow the existence of a densely-defined inverse, but an operator 
providing the inverse as the classical limit can nevertheless be defined ||. When one 
departs from the classical regime, however, quantum corrections arise whose form can 
be computed in some models and which can be parametrized sufficiently generally 
for phenomenological investigations. In the presence of these corrections, provided 
they are small, consistent gauge-invariant cosmological perturbation equations have 
been determined 0. 

Gauge invariance in general relativity and the candidates for its quantization 
is intimately related to spacetime structure because gauge transformations include 
changes of coordinates. Accordingly, quantum corrections to gauge-invariant pertur- 
bation equations show how fundamental quantum spacetime effects such as discrete 
geometry are reflected, in physical terms, by implications for cosmological observ- 
ables. Cosmology then provides an intriguing link between the fundamental physics 
of spacetime and phenomenology, as it will be explored in this paper. 

Looking at the issue from the phenomenological side, linearization of dynamical 
equations in the presence of inhomogeneous perturbations is one of the most exten- 
sively studied problems in cosmology. With recent progress in LQC, the context in 
which such questions can be addressed has been extended to a new class of quantum- 
gravity effects with the aim of making a prediction for early-universe spectra and, 
hopefully, constraining the theory. The perturbed equations contain quantum cor- 
rection functions and are augmented by counterterms which guarantee cancellation 
of anomalies in the effective constraint algebra [0, H 3 . These contributions provide 
insights into the spacetime structure, modify the equations of motion and, eventually, 
imply characteristic signatures for physical observables. 
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The perturbed equations of motion for scalar, vector, and tensor modes have 
been computed, respectively, in |6|, [7| ||, while the tensor spectrum and index have 



been found and explored in |10|, |TT[| . To close the set of cosmological consistency 
relations, it remains to compute the scalar spectrum and its derived observables. 
Such is the first goal of this paper: We shall find the Mukhanov equation of scalar 
perturbations for the first time and solve it. With exactly the same procedure, we 
will rederive the solution of tensor modes with zero effort. The scalar and tensor 
spectra and spectral indices will be obtained together with a consistency relation 
between the tensor-to-scalar ratio and the tensor index. In doing so, we will derive a 
conservation law for the curvature perturbation, extending the well-known classical 
result to the quantum-corrected equations, and discuss its implications for quantum 
spacetime structure. 

All cosmological observables, including index running and higher-order quanti- 
ties, are linearly corrected by quantum terms Sp\ whose magnitude we cannot predict 
yet because we have no control over the details of the underlying full quantum theory. 
The presence of extra parameters would seem to at best make it possible to place 
upper bounds on the quantum corrections 5p\ for a given inflationary potential. For 
instance, one can arrange to have large-enough quantum corrections so that the scalar 
running be sizable. On the other hand, we can naturally envisage a situation where 
5p\ is much smaller than the slow-roll parameters, and therefore completely negligible. 
Until better control over the theory is achieved, the phenomenological consequences 
of inverse-volume (or other) corrections will remain unclear but here, as our second 
goal, we highlight the following issue: Different parametrization schemes will lead to 
different background solutions and predictions for the size of Sp\. In particular, the 
usual minisuperspace parametrization of FRW loop quantum cosmology seems to 
be incompatible with anomaly cancellation in inhomogeneous LQC, as well as with 
the simplest power-law solutions. Conversely, the lattice refinement parametrization 
overcomes all these problems, can predict much larger quantum corrections, but it 
might indicate a problem related to superluminal propagation of the perturbations. 
We wish to stress all these features and the importance of further investigating the 
lattice refinement parametrization, which requires input from the full theory. 

Since anomaly cancellation so far has been shown to occur only in the quasi- 
classical regime where inverse- volume corrections and counterterms are small, we 
shall concentrate on this case. Therefore it is not possible to draw comparisons with 



holonomy-correction results |10|, [12], |13[ or with previous works which partially fixed 
the gauge, considered test fields and neglected metric back-reaction, since they were 
devoted to the superinflationary regime of a near-Planckian epoch [[14], [15], [16], [L7|, [18], 



19]. 



In section |] we review the LQC background equations of motion for a flat FRW 
model and a scalar field, and the parametrizations arising in minisuperspace quanti- 
zation and the lattice refinement phenomenological approach. Background solutions 
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with de Sitter and exact or perturbed power-law expansion are found in section |[ 
Section || is completely devoted to scalar perturbations. After a review and some 
important updates of the results of 0, we show that the comoving curvature per- 
turbation is conserved at large scales, as in the classical case (section |5.3| ). The 
Mukhanov equation for a scalar perturbation variable is then found and solved in 
section |5.4|, while the set of linear cosmological observables (power spectrum, index 
and index running) is derived in section p.(5\. The set of observables is completed 



in section |], where the same analysis of the previous section is applied verbatim to 
tensor perturbations. A discussion of the main achievements of this paper and future 
directions can be found in section [f[ 

2. Spacetime structure and phenomenology 

Loop quantum gravity has provided results that show the discreteness of spatial 
quantum geometry: geometrical operators such as those for area and volume have 
discrete spectra |2D|, PT| . Taken by themselves, these features are not observable be- 



cause the corresponding objects, the areas of spatial surfaces or the volumes of spatial 
regions, are not gauge invariant. However, these mathematical properties, derived 
from the underlying principle of background independence on which the theory is 
built, indicate new features of the quantum representation with important effects for 
the dynamics. The volume operator, for instance, enters matter Hamiltonians 
and the gravitational Hamiltonian constraint operator |22| , and thus influences their 
properties at the quantum level. In this way, once the dynamical equations are suffi- 
ciently well understood, physical observables are affected and a potential comparison 
with observations is made possible. 

The dynamics of loop quantum gravity amounts to that of a coupled, interact- 
ing many-body problem in which the elementary constituents are the fundamental 
building blocks of space. Such equations governing the dynamics are difficult to solve 
exactly, but several crucial effects visible in them are generic and characteristic; they 
provide the basis for phenomenological evaluations. There are two main effects: (i) 
inverse-volume corrections and (ii) holonomy corrections due to the fact that the 
background-independent quantization used in LQG allows the representation only 
of exponentiated curvature components, gathered as the holonomies of connection 
variables. Inverse-volume corrections are currently under much better control in 
cosmological perturbation equations, and so they will be our main focus here. 

To see implications of those effects in the dynamics, the many-body Hamiltonians 
must be analyzed. There is by now a systematic procedure to do so, based on effective 



canonical equations to describe semi-classical dynamics |23], |24j. Effective equations 
in this context, analogous to low-energy effective actions for expansions around the 
ground state of anharmonic systems, are obtained from expectation values of the 
Hamiltonian (constraint) operators in a generic class of semi-classical states. So far, 
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these equations have been computed and analyzed only in isotropic cosmological 
models, in which a solvable system analogous to the harmonic oscillator is available 



25]]. The form of the equations, however, is more general and can be used also in 
the presence of inhomogeneous cosmological perturbations. To avoid bias, one only 
has to ensure that correction functions are parametrized sufficiently generally, since 
the control over the theory is not yet strong enough to provide unique predictions 
for them. 

These equations, including parametrized quantum-gravity corrections, allow im- 
portant conclusions at the fundamental and phenomenological levels by a clear line 
of arguments. First, the structure of spatial geometry changes according to LQG, 
affecting the form of the constraints generating gauge transformations. Secondly, in 
the presence of quantum corrections the algebra under Poisson brackets (or commu- 
tators when quantized) of the constraints as gauge generators is modified, amounting 
to a different realization of the classical transformations of spacetime and the under- 
lying gauge behaviour. This abstract feature has several consequences. For instance, 
quantum-geometry corrections cannot be implemented by any higher-curvature ef- 
fective action because those corrections would not change the classical constraint 
algebra and the underlying notion of gauge and covariance. An effective action that 
can include all effects of LQG must be of a more general form, for instance allow- 
ing for non-commutative geometry. Moreover, gauge transformations belonging to 
a deformed algebra no longer correspond to ordinary coordinate transformations on 
a manifold. Thus, effective line elements may be questionable because the transfor- 
mations of metric components and coordinate differentials collected in ds 2 no longer 
match to make the line element invariant. In such a context, physical information can 
be gained only from gauge- invariant variables that take into account the new gauge 
structure. With these additional effects, there is a chance that quantum gravity cor- 
rections may be stronger than usually expected, for instance, from naive arguments 
based on the size of higher-curvature corrections. 

One possible new phenomenon is of particular interest. When the constraint 
algebra is deformed, the Bianchi identity or the conservation equation for stress- 
energy is modified. There are still analogous identities if the deformation is consistent 
and anomaly-free, but they may refer to different quantities than in the classical 
case. Then, the curvature perturbation is no longer guaranteed to be conserved, a 
possibility which has already been raised R Ef|. If the curvature perturbation is 



no longer conserved, on the other hand, magnification effects for modes outside the 
Hubble radius can be expected. Even though deviations from conservation given 
by quantum gravity were small at any given (sufficiently late) time, expected to be 
determined by the tiny ratio of the Planck length by the Hubble radius, the lever arm 
of non-conservation during long times between horizon exit during inflation and re- 
entry might magnify those tiny effects. In this way, a tight link is obtained between 
fundamental spacetime structure and cosmological phenomenology. 
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Here we will provide results in both directions, fundamental physics as well as 
phenomenology. In particular, we will demonstrate that a subtle cancellation in 
the anomaly-free correction functions of LQC does make the curvature perturbation 
conserved in spite of the non-trivial deformation of the constraint algebra, a feature 
which has not been noticed before. As a consequence, effective linear perturbations 
of Friedmann-Robertson- Walker geometries can be meaningfully constructed even 
in the presence of a deformed gauge structure. Corrections to standard perturbation 
equations then follow naturally. 



3. Background equations and parametrizations 



To begin, we write down the LQC effective equations of motion for an FRW back- 
ground ds 2 = a 2 (r)(— dr 2 + dxMxj) in conformal time r (see, e.g., the review in 
11] for a detailed derivation of these results). We shall ignore holonomy corrections, 



which have not yet been considered in the perturbed dynamics. Notice, however, that 
these contributions in some parametrizations dominate over inverse-volume correc- 



tions as far as tensor modes are concerned [11], [L3|; we will further discuss this point 
in the final section. For a scalar field tp with potential V, the effective Friedmann 
and Klein-Gordon equations read 

j2 



8nG 



and 



<p" + ZH[l- 



2u 



+ pV(<p) 



d In v 
dlnp 



<f' + vpV !(p = 



(3.1) 



(3.2) 



where G is the gravitational constant, T-L = a' /a is the Hubble parameter, primes 
denote derivatives with respect to conformal time, and p = a 2 (in comoving volume 
units) is the triad variable in minisuperspace. (Triad variables can take both signs 
depending on the orientation of space. Here we assume p to be positive without loss 
of generality for effective equations.) From equations ( |3.1| ) and ( |3.2| ) one obtains the 
Raychaudhuri equation 



In these equations, 
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and a and pp\ = a|j are constant. 1 We will often need to switch from p to conformal 
derivatives via 

' = (3.7) 
amp 

and the formulae 

5' Pl = -am PU 5^ = aH 2 (a -l + e)5 PU (3.8) 

with e = 1 — %' /% introduced as a slow-roll parameter below. 

Functional forms of a(p) and v(jp) are fully computable in general form from op- 
erators in exactly isotropic models |27| and for regular lattice states in the presence 
of inhomogeneities [p8|] , with parametrizations of quantization ambiguities affecting 
the values of ao, v$ and a [^, |30||. However, only the expanded forms ( |3.4|) and (|3.5| ) 
are needed in the perturbative regime considered here. From these explicit calcula- 
tions of inverse- volume operators and their spectra, one can derive further properties 
characteristic of loop quantum gravity. In particular, correction functions imple- 
menting inverse-volume corrections, when evaluated at large values of the densitized 
triad or the scale factor in a nearly isotropic geometry, approach the classical value 
always from above. This consequence, which is a robust feature under quantization 
ambiguities and will turn out to be important later, implies that the coefficients ao 
and pq introduced in the parametrizations used here must be positive. 

3.1 Parametrizations 

The above equations are derived in a minisuperspace Hamiltonian formalism where 
the super-Hamiltonian (the only non-trivial constraint on homogeneous backgrounds) 
is first symmetry-reduced, and then quantized with LQG techniques. The resulting 
equations then constitute partial effective equations, which means that they cap- 
ture the behaviour of expectation values of observables in semi-classical states, but 
without taking all quantum corrections into account. In particular, quantum back- 
reaction by fluctuations and the holonomy corrections of loop quantum gravity are 
not included at the present stage. The LQG techniques consist in a choice of canon- 
ical variables and operator ordering which make the final result quite different with 
respect to the traditional minisuperspace Wheeler-DeWitt quantum cosmology. At 
the semi-classical level, the main difference is the presence of correction functions 
and (|3.5|). The constants ao, vo an d cr will enter the cosmological observables 



and it is important to set their value range beforehand. This range strongly depends 
on the physical interpretation of the model. We can identify two views on the issue, 
one purely homogeneous and isotropic and the other associated with the natural 
presence of inhomogeneities. 



1 We put a subscript 'PI' in the definition (3J3) in order to avoid confusion with perturbations 



such as dip. However, the equations below do not rely on any particular characteristic scale dpi, 
which may differ from the Planck length. 
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3.1.1 Minisuperspace parametrization 

On an ideal FRW background, open and flat universes have infinite spatial volume 
and the super-Hamiltonian constraint is formally ill defined because it entails a di- 
vergent integration of a spatially constant quantity over a comoving spatial slice 

J d 3 x = +00 . 

To make the integral finite, it is customary to define the constraint on a freely chosen 
finite region of size V = a 3 Vo, where Vo is the corresponding comoving volume: 

d 3 a: ->■ / d 3 x = V < +00 . 

S ^S(Vo) 

The volume appears in the correction function (|3.6| ) as 5-p\ ~ a _cr ~ V~ a/3 . To make 
5pi adimensional, one can use the Planck length £p\ to write 

'PI 



5pi ~ {fj a- CT . (3.9) 

To specify the coefficient further, one sometimes introduces the area gap Api = 
2y/3n'j£p l , where 7 is the Barbero-Immirzi parameter. 2 A detailed calculation then 
shows that the constant coefficients ao and Uq are fTl] 



(3q-o)(6q-o) /ApA 2 o(2- I) / A P1 \ 2 

« = ^ (— j • = -jj- (— J • P-IO) 

They depend on two sets of ambiguities, one (1/2 < / < 1 and 1/3 < q < 2/3) 
related to different ways of quantizing the classical Hamiltonian 3 and the other (a) 
depending on which geometrical minisuperspace variable has an equidistant stepsize 
in the dynamics: in terms of the triad variable p, p a ^ 2 is equidistant if inverse-volume 
corrections with exponent a in Api appear. More physically, this parameter is related 
to how the number of plaquettes of an underlying discrete state changes with respect 
to the volume as the universe expands. The latter is a phenomenological prescription 
for the area of holonomy plaquettes, but ideally it should be an input from the full 
theory |34 . 



In the minisuperspace context, a natural choice of these parameters is 

a = 6, 1 = 1, q = \, (3.11) 



2 The use of twice this value according to recent findings may be better justified |3T| but the 
resulting change in the values of ao and vq is not relevant for what follows. 

3 The different interval for q with respect to the one given in |llf| stems from the same argument 
in the full theory which constrains the range of I ^2|. We set the natural value of q to be 1/2 



rather than 1 [B3| in equation (3.11) 
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so that, assuming pp\ = Api [33| 



, one has 



« = ^~0.04, u = |»0.14. (3.12) 

Notice that o = 6, motivated by holonomy corrections not becoming large at small 
curvature, corresponds to the so-called 'improved quantization scheme' p3| , a name 



which applies also to values of I and q different from equation (|3.11| ). A specific 



example provided in p3| , for instance, is q = 1 such that a = 0. Note, however, that 
the exact value a = is obtained in this case only due to a spurious cancellation 
in isotropic settings; it seems to suggest negligible inverse-volume corrections but 
is unreliable compared with more general derivations. For phenomenology at the 
current level of precision, the most significant parameter is a, which is not as much 
affected by different choices of the minisuperspace scheme. 

Since 5p\ is Vo-dependent, inverse- volume corrections cannot strictly be made 
sense of in a pure minisuperspace treatment. Inverse- volume corrections, as used 
here, have never been derived fully consistently in this context owing to the In- 
dependence. In particular, as discussed in [^5], [36[, while the improved scheme does 
take into account refinement for holonomies in an ad-hoc manner, it ignores these 
effects for inverse- volume corrections. This failure to represent inverse- volume effects, 
which are crucial for well-defined Hamiltonians in loop quantum gravity, presents 
a serious limitation of pure minisuperspace models which can be overcome only by 
bringing in further ingredients to take into account the behaviour of inhomogeneities, 
as indicated in what follows. Precise derivations become more complicated in this 
situation, but by a combination with input from phenomenology one can obtain 
valuable restrictions on the possible realizations. 

We wish to make a further comment on this issue. In the context of inflation, 
one has a set of observables given by the anisotropy spectra and their derivatives; all 
these quantities contain parameters of the inverse volume corrections. Since they are 
all evaluated at horizon crossing, the comoving fiducial volume therein (implicitly 
conceived as greater than the causal region ~ "H~ 3 ) can be naturally set to be the 
Hubble volume. If one maintained the minisuperspace parametrization also in the 
presence of perturbations, the conclusion would be that the fiducial volume problem 
is less severe than expected. However, the minisuperspace setting pertains only to 
exactly isotropic models, and this solution of the problem is at best incomplete. To 
get a clearer picture, we should include inhomogeneities already at the fundamental 
level. The following argument shows how to do so qualitatively. 

3.1.2 Lattice refinement parametrization 

The chosen volume V is a purely mathematical object which should not appear in 
physical observables, but it does appear in equation (|3.9|). Since 5p\ will enter the 
observables, we might face a problem. To make the situation better behaved, we 
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introduce generic effects of inhomogeneities. One example for doing so is the lattice 



parametrization discussed in [11| which, as one implication, extends the range of the 



parameters. Then, one has a large range of a, 

4<a<6, (3.13) 

the value of the improved minisuperspace quantization scheme being included. When 

Pp\ = Api, for instance, varying over the range of a (equation ( p.!3| )), I and q (as 
specified below equation (|3.10|) ), 

0<a <^~0.06, (3.14) 

0.07^ % < i/ < | ~0.17. (3.15) 



Instead of repeating the arguments leading to equation ( |3.13| ), we consider an 



alternative lattice parametrization where one uses the 'patch' volume of an underlying 
discrete state in correction functions, rather than the much larger volume V fl36fl . 
In this parametrization, motivated by key aspects of discrete spacetime dynamics, 
the ranges of parameters change more significantly than in ( j3.13j ), with important 
consequences for phenomenology. Now, corrections refer to the patch size 

v = ^, (3.16) 

with Af, the main input from quantum gravity, the number of discrete patches in V. 
By construction, v is independent of the size of the region, since both Vo and Af scale 
in the same way when the size of the region is changed. Physical predictions should 
not feature the region one chooses unless one is specifically asking region- dependent 
questions (such as: What is the number of vertices in a given volume?). 

Given the available parameters and their dimensions, the leading-order quantum 
correction in a and v is then of the form 



*P1 1 I /?3 • /V 



>!•! 1 — j = \ t P1 — 



(3.17) 



where m > is a positive integer parameter. It determines the power by which p 
appears in leading corrections of an expansion of inverse- volume correction functions. 
Primarily, the correction functions a and v, and thus dpi as well, depend on flux 
values, corresponding to p for the isotropic background. Since p changes sign under 
orientation reversal but the operators are parity invariant, only even powers of p 
can appear, giving m = 4 as the smallest value. At this stage of development of 
the full theory, it is not entirely clear that general correction functions depend only 
on the fluxes (rather than, e.g., also on the eigenvalues of more complicated volume 
operators; for properties of their spectra see |J7], |3B|). Therefore we set m > 4. 
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A time-dependent J\f(t) corresponds to the dynamical 'lattice-refinement' be- 
haviour []34| . For some stretches of time, one can choose to use the scale factor a as 
the time variable and represent J\f(a) as a power law 

A/" = A/" a- 6x , (3.18) 

where Ao is some (coordinate and Vo-dependent) parameter and the power x de- 
scribes different qualitative behaviours of changing lattices. Overall, we have 

/ r \ — 

<y pl = (4^) 3 a-^ m . (3.19) 



This equation cannot be obtained in a pure minisuperspace setting (e.g., fill ) where 
only one parameter a enters as in (|3.9| ). The presence of an extra parameter, com- 
pared to the minisuperspace parametrization, may appear as a disadvantage, but 
we will see later that it is required for being able to match with phenomenology; 
otherwise the model would be ruled out. 

In the lattice-refinement derivation, the parameter a plays two roles, one as a 
dynamical geometric quantity and the other as internal time. While writing down 
the semi-classical Hamiltonian with inverse-volume (and holonomy) corrections, one 
is at a non-dynamical quantum-geometric level. Then, internal time is taken at a 
fixed value but the geometry still varies on the whole phase space. In this setting, 
we must keep M fixed while formulating the constraint as a composite operator. The 
net result is the Hamiltonian constraint operator of the basic formulation of loop 
quantum cosmology |39, |4(| not taking into account any refinement, corresponding 



to x = and m = a. However, when one solves the constraint or uses it for effective 
equations, one has to bring in the dynamical nature of M from an underlying full 
state. This is the motivation for promoting M to a time- dependent quantity, a 
step which captures operator as well as state properties of the effective dynamics. 
Its parametrization as a power law of the scale factor is simply a way to encode 
the qualitative (yet robust, see below) phenomenology of the theory. The general 
viewpoint is similar to mean-field approximations which model effects of underlying 
degrees of freedom by a single, physically motivated function. 

Comparing with the earlier minisuperspace parametrization, equation Q3.19] ) 
gives a = (2x + l)m as far as the a-dependence is concerned. The number of 
vertices TV must not decrease with the volume, so x < 0; it is constant for x = 0. 
Also, v = V/ftf ~ a 3 ( 1+2:r ) is the elementary geometry as determined by the state; 
in a discrete geometrical setting, this quantity has a lower non-zero bound which 
requires — 1/2 < x < 0. In particular, for x = —1/2 we have a constant patch 
volume, corresponding to what is assumed in the improved minisuperspace quanti- 
zation scheme p3|. In contrast with the minisuperspace parametrization (|3.11|), in 



the effective parametrization of equation Q3.19 ) we have a = in this case. Thus 
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even for x = — 1/2 is the parametrization new and different from the minisuperspace 
representation, overcoming the problem of representing inverse-volume effects in a 
pure minisuperspace treatment. 

To summarize the general lattice-refinement scheme, o = (2x + l)m is a time- 
independent or slowly changing parameter 4 given by the reduction from the full 
theory and with range 

<7>0. (3.20) 

Assuming this range will be of utmost importance for justifying the validity of the 
cosmo logical perturbation spectra. In fact, we shall find that a must be small in order 
for the spectra to be almost scale invariant, a range of values that cannot occur for 
either the minisuperspace parametrization or the first lattice parametrization ( |3.13| ). 
Responsible for the better matching is the new parameter x, while m alone (or x — 0) 
would give a range similar to the minisuperspace parametrization. 

Note that, in principle, a may be different in a and v for an inhomogeneous 
model. However, here we assume that the background equations fl3.4Q and (|3.5Q with 
the same b-p\ are valid also in the perturbed case. Not only is this choice natural 
whenever background quantities are considered, but it is also crucial for several 
simplifications to follow. 

Before moving on, a remark is in order. The patches of volume v find a most 
natural classical analog in inhomogeneous cosmologies, in particular within the sepa- 
rate universe picture jll]] . For quantum corrections, the regions of size v are provided 
by an underlying discrete state and thus correspond to quantum degrees of freedom 
absent classically. However, the discrete nature of the state implies that inhomo- 
geneities are unavoidable and no perfectly homogeneous geometry can exist. Given 
these inhomogeneities and their scale provided by the state, one can reinterpret them 
in a classical context, making use of the separate universe picture. There, the volume 
V can be regarded as a region of the universe where inhomogeneities are non-zero 
but small. This region is coarse grained into smaller regions of volume v, each cen- 
tered at some point x, wherein the universe is FRW and described by a 'local' scale 
factor a(t, x) = a x (t). The difference between scale factors separated by the typical 
perturbation wavelength |x' — x| ~ A <C V 1 ' 3 defines a spatial gradient interpreted 
as a metric perturbation. In a perfectly homogeneous context, v ~ V and there is 
no sensible notion of cell subdivision of V; this is tantamount to stating that only 
the fiducial volume will enter the quantum corrections and the observables, M = AV 
On the other hand, in an inhomogeneous universe the quantity v carries a time de- 
pendence which, in turn, translates into a momentum dependence. The details of 
the cell subdivision (number of cells per unit volume) are intimately related to the 

4 The creation or subdivision of new cells in a discrete state depends on the spatial geometry and 
can thus be considered as changing more slowly than other processes in an expanding universe. On 
large time scales, the parameter a may change, distinguishing different microscopic epochs in the 
history of the universe. 
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structure of the small perturbations and their spectrum. Thus, lattice refinement is 
better suitable in the cosmological perturbation analysis. As long as perturbations 
are linear and almost scale invariant, the size of volume within which the study is 
conducted is totally irrelevant. 

3.1.3 Parameter estimates 

Since the lattice refinement picture is phenomenological, presently we are unable to 
determine the quantum correction Sp\ from first principles and, at this stage, the 
latter is regarded as a free parameter which can be constrained by experiments. 
In the minisuperspace and lattice parametrization (|3.13|) , on the other hand, there 



may be an argument which estimates the magnitude of dpi heuristically. In fact, in 
those parametrizations 5p\ = (5o/V) a ^ 3 , where So is some constant volume and V is 
the fiducial volume; for lack of better knowledge, one often assumes <5o ~ £ Pl . All 
inflationary observables are evaluated at horizon crossing, so the volume V is very 
naturally fixed by the size of the Hubble horizon at that moment (denoted with a *) 

V~#~ 3 , (3.21) 

where H = a/ a = H/a. (This equation is invariant under isotropic rescalings of 
the coordinates. In terms of dpi introduced via 6p\ = (api/a) a in the minisuperspace 
parametrization, we may write ap\/a = £p\H*.) The point here is that so far V 
has been arbitrary, the only requirement being that it contains the Hubble region 
at any given moment. Provided 5q is fixed a priori, this equation fixes V once and 
for all because one is not at the liberty of changing the numerical factor in ( |3.21| ), 
which is 0(1). Slightly different definitions of the Hubble horizon differ only for 0(1) 
coefficients, which do not affect the discussion qualitatively (on the other hand, 0(10) 
or O(O.l) coefficients are unacceptable because the observables here are defined at, 
not before or after, horizon crossing). 

To estimate 5p\ during inflation, we could take the grand-unification scale ~ 
10 14 10 17 GeV. As 5p\ = (4n#*) ff and 4 < a < 6, one has the upper bound 
5p\ < O(10 -8 ). Typical prefactors in observable quantities may even carry an extra 
O(10 _1 ) suppression, as we will find, so none of the inverse volume corrections with 
these choices are observable, even in the scalar running, as the slow-roll parameters 
are O(10~ 2 ). In comparison, the inflationary tensor index gets an extra contribution 
<5hoi oc (£p\H*) 2 < O(10 -4 ) from holonomy corrections [fL3| , which dominate over Sp\. 



Holonomy corrections in the perturbed scalar sector has never been computed, but 
we expect a similar hierarchy of scales. 

Unfortunately, these estimates rely on a particular choice for 5q and the argu- 
ment cannot be regarded as robust. The size of this dimensionful constant strongly 
depends on the underlying theory, and a change in magnitude of 5 would modify the 
above results. Moreover, we will see that the minisuperspace parametrizations are 
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not compatible with requirements on background solutions combined with anomaly 
freedom. The estimate of 5p\ could thus at best be used as an external input for the 
lattice-refinement parametrization. Then, once the scale of inverse-volume correc- 
tions is fixed, a further consistency condition must be satisfied because inverse-triad 
corrections and holonomy corrections are determined by the same parameter that 
specifies the underlying discreteness scale []5B| . To evaluate this condition, we use 



the alternative form of 8p\ = (£ Pl /v) m ^ 3 = ipjv 4 ^ 3 for m = 4. The underlying dis- 
creteness scale, as a distance parameter, is then the linear dimension of patches, 
L = v 1 / 3 = ^pi/5p{ 4 = fpi/^pi-f/*)' 7 / 4 . On the other hand, the strength of holonomy 
corrections can be expressed in terms of the critical density p cr i t = 3/(87rG7 2 L 2 ); 5 
holonomy corrections are weak when the matter density satisfies p <C Pcrit- With L 
as assumed here, 

Pcrit = „ o ~ P* \ — , (6.22) 



87rG7 2 \Ppi 

with the matter density p* at the time of horizon crossing. Thus, for p crit 3> p* to 
ensure that holonomy corrections do not significantly alter the classical behaviour at 
horizon crossing, we must require a < 4. Then, 5p\ becomes larger than estimated 
above. For a critical density of Planckian size, which is often desired so as to have 
strong quantum-gravity corrections only in the Planckian regime, a must be close to 
zero. For such values, 5p\ ~ 0(1) with the above estimate of 5 , clearly dominating 
holonomy corrections at horizon crossing. (For a = 2, the critical density is the 
geometric mean p crit ~ y/pp\P* and we have 8p\ ~ <5hoi-) 

3.2 Slow-roll parameters 

For later convenience, we define the first three slow-roll parameters as 

(3.23) 
(3.24) 

+ e + ri-l, (3.25) 





w 


1 - 


H 2 ' 


1 - 


</ 




Hip' 


1 


K 


w 


w 



which coincide with the standard definitions in synchronous time 

H 2 ' '~ Rip' ^ ~ H 2 \<p 



5 As an explicit calculation shows, the true critical density is actually ap cr it pl| . Setting q = 1 
and a — 6 in the lattice parametrization one gets a = 1, as in j33j. In general a does appear, but 
the arguments in this discussion are qualitative and we can ignore this issue. 
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The parameter e will be especially important later on and we can rewrite it as 

v-' / ldlnzA d In a 



1 



3 dlnpy dlnp 



12 



AttG 



f_ 



i + 



a + u 



a 



6 



1 



8 



pi 



, cra 



(3.26) 



using the Raychaudhuri equation ( |3.3|) in the first step. Notice that the symbol = 
in the last line of equation ( |3.26 ) implicitly hides the 0(5pi) truncation. This note 
of caution applies to any of the equations below, where 0(5 P1 ) terms are dropped 
as required for self-consistency of perturbed equations. In contrast, the slow- roll 
approximation will always be invoked explicitly and indicated with the symbol ~. 
The derivatives of e and r\ are 



2H[e + 



-2H 



d In a 
dlnp 

d 2 In a 



7] + 



d In a d In v 1 d 2 In v 



dlnp d\np 3d hip 2 



I 1 d In v 
3 dlnp 



rj 

where 



din p 2 
2Ut{e - rj) - 



a'Ktb 



pi 



a - + 2e - rj 



a 

- - 1 

6 



(3.27) 
(3.28) 

(3.29) 



While in standard inflation e is almost constant whenever it is small (since the 
classical part of e' is quadratic in the parameters), depending on the size of the 
quantum correction ae5py the quantity e' could be of the same order as e. However, 
we expect 5pi to be small in the typical setting. 

For a given background a(r) and <p{t), the slow-roll parameters are functionally 
identical to the classical case. Clearly, the potential required to give rise to such an 
evolution is different, as one can see also from the (later useful) relations 



Up 



V,tptp — 



up 
U 2 

Up 

1 



r?-3 + 2 
3(e + 77) 



dlnu 
dlnp 

au 5 Pl ) 



(3.30) 



V 



e + 2(3 - e 



, d In u H 2 In u 
amp amp z 



dlnz/ 
dlnp 



up 



(-m 2 + W 2 aii tp dpi) 



where 



^V + ( 2 -3(£ + i))] 



fi v = u (a - 3 + e + 2r]) 



(3.31) 

(3.32) 
(3.33) 
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4. Background solutions 



Let be a set of generic scalar variables and let us write the background equations of 
motion, as well as the soon-to-be-found Mukhanov equation for the scalar perturba- 
tion, as 0[(p} = 0, where O is a (possibly non-linear) differential operator. One can 
drop quantum terms of order higher than 8p\ and split each variable into a classical 



part 4> c and a quantum correction q 5pi |11 



> q opi 



so that each equation becomes 



a[0 c ] + {a[0 q ] + aj0c]}5pi = o 



(4.1) 



(4.2) 



Requiring that the classical and quantum part vanish separately (a condition which 
defines what is meant by (f> c ) yields two equations: 



a 



o, 



a 



0. 



(4.3) 



This splitting strongly resembles the one into coarse- and fine-grained perturbations 
in stochastic inflation H2, |3, 31], E5J; in fact, for a Klein-Gordon scalar the sec- 
ond equation ( (4.3| ) is nothing but a 'complementary' Langevin-type equation for a 
quantum variable with a noise term sourced by the classical part. 
For example, consider the scalar field and scale factor profiles 



ip = ip c + ip q 5pi , a = a c + a q 5pi . 
The Hubble parameter can be written as 



(4.4) 



H — T-L c + HqSpi , 



- (1 + <x)H c 



(4.5) 



Also, the scalar potential V is expanded in a Taylor series around (p c , 

V(<f) = V(if c ) + V^(<£ c )<£ q 5pi = V c + V q 5pi , 

V lV (tp) = V tV (<p c ) + V w (Lp c )Lp q 5pi = V vc + V tVq 6pi . 



(4.6) 
(4.7) 



Plugging these expressions into the Friedmann and Klein-Gordon equations (3.1) 
and (p.2|), we obtain a pair of classical equations, 



C c I ' 



= rt + ZH c <p' B + alV tVC , 



(4.8) 
(4.9) 
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plus another pair of relations involving the correction functions (f q and a q : 



4ttG 



"0-^0 ,2 



+ a c > K + Kj) + 2a q V c 



= <p'Z+ 2H C (1 - a>' q + aH 2 c (a + e c - 3)<p q + a 2 c (v V vc + V >vq ) 
+2a q V vc + cp' c {2n q + av Q U c ) . 



,(4.10) 



(4.11) 



At this point we look for special background solutions with exactly constant 
slow-roll parameter e, i.e., with a scale factor expanding as a power-law: 



a = a c = r 



n< -1 



(4.12) 



where r < and the limit n ~ — 1 corresponds to de Sitter spacetime. By this 
definition, the quantum corrections a q and "H q vanish identically. At the classical 
level, one gets the power-law solution [f46| 



(4.13) 
(4.14) 



n(n + 1) 

ip c = Lf In \t\ = ±\l — : — - — In |rl 



y = Y Q-Hn+Vp/va 

for which V c = V(ip c ) = Vo\t\~ 2ti ~ 2 and 



A-kG 

n ( 2n ~ l ) -2(n+l)vM) 

8ttG 



e = Vc = £ c = 1 



n 



(4.15) 



Let us see if there exist solutions of the form ip q = ip q0 \r\ b , distinguishing two cases: 
• If n ^ -1, 

V n = -2(n + l)V n ^\r\ b - 2n - 2 , (4.16) 



(j -2{n + l)V^\r\ b - 2n - 2 



for Vq as defined in (|4.6|) . Equations ( |4.10[) and (|4.11| ) become 



= WqO 



-an)-2(n + l)-| 



= <^q0 



«0 - Vq 
~ h «0 — 

2 y? 2 



:^0 



T 



(4.17) 



6(6 - 1) + 26n(l - a) + an(an + 1 - 2n) + 4(n + l) 2 ^ 



T 



1 6—2 



an — 2(n + 1)— | 



(4.18) 



It turns out that, if (p q0 ^ 0, the solution requires 6 = and the equalities 



3n 

1 +71 



a - Vq 



2(an + 2n- V 



Vq a Q 

e 



1-e 



2(3 + a -e) 



(4.19) 



a 



= vq [6(3 - e) - <r(3 + o- - e)] - 3a 2(3 - e) - -(3 - a - e) . (4.20) 
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The potential then reads 



V 



l_2( n + l)^ P1 



V, 



while the second and first slow-roll parameters are 



2 ^qO r 

T) = n c -a n dpi 



f = £-<7>n + l + n)^ P1 . 



An exact solution is 



a 



0, 



a = 



^=(l-a 5pi)K, 



(4.21) 

(4.22) 
(4.23) 

(4.24) 



while for a > 0(1) and small e the expression ( |4.20|) is satisfied if a and v§ 
obey 



6 



7< 3 



0. 



(4.25) 



[18 - a(3 + o-)] 

For < a < 3, the solution has «o, z/ > if e > cr(3 — cr)/6 (e.g., e > 1/3 if 
a = 1 or a = 2). Because of the lower bound on e, this solution prefers the 
limiting values a ~ 0, a ~ 3 if extreme slow roll is to be realized. 

When a > 3, a and z/ have opposite sign and |i/ | ^ |«o|- This case is excluded 
by the above considerations on inverse-volume operators, which require «o and 
b>o to be both positive. 

If y? q o = 0, equations ( |4.17| ) and ( |4.18| ) are solved for 



a 



V 



(4.26) 



This solution has z/ ^> «o an d 2 < a < 3. For a given background, e is constant 
per ( |4.15|) ; a relation to the constant a , an d cr may thus be acceptable. 
However, the required tuning of general quantization parameters to background 
parameters makes this solution very special. 

By construction and consistently, all these power-law solutions obey equation 
( |3.26|) . Their qualitative features are summarized in table ||. 

The last exact power-law case we consider is de Sitter, n = — 1, whose classical 
solution is tp c = const. There we cannot use equation (|4.16|) because the 
potential (jUg ) was derived using ( |4.13| ) which is ill-defined for n — — 1. From 
equation ( |4.10| ) we simply obtain V q = — a V c . Since V q = V v (<p c )(p q , this 
implies that <p q is constant. But 5p\(a) is not constant, so that ip is not constant 
and V(ip) can be reconstructed from the evolution. Combining ( |4.1|) and (fO|), 
we find 



V((p) = K(l - u 5pi) = V C + V jtpc (<p - <p c ) 



(4.27) 
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a 


a ,v 


<Pc(t) + V^qO^Pl 


a = 

< a < 3 
a > 3 


a = u 
oto, Vq > if e > cr ( 3 ~°") 
sgn(a ) = -sgn(fo), \u \ » |a<o 


<Pc(t) 


2<a<3 


«o = f ^0 < ^0 



Table 1: Inflationary power-law solutions a = |r| n with 0<e = l + l/ra<l and 
exponential potential. ip c (t) and ip q o are given by equations (|4.13| ) and Q4.19| ), respectively. 



Equation (|4.11| ) becomes a (a 



3)T VqO + O VqO 



0, where m 2 = V %iptp {(p c 



The reconstructed potential above is linear in <p, so m 2 
implies that either o = or o = 3. 



0. If <^ q o ^ 0, this 



In the next section we will find that not all of these solutions will be compatible 
with a certain consistency relation on quantum counterterms. We have not shown 
that the above solutions are attractors in configuration space, a necessary condition 
for adopting them as valid backgrounds. In the next section we will assume this 
is the case, since in the quasi-classical regime the dynamics is very close to general 
relativity. Anyway, the structure of perturbation equations and observables does not 
change qualitatively if one expands about a more general quasi-de Sitter solution. 

Also, setting a = a c as in equation Q4.12D will not result in any loss of generality 
in solving the Mukhanov equation. Since the coefficients in the second equation (|4.3| ) 
will depend only on T-t c , the structure of the quantum corrections in the solution will 
be always the same, regardless of "H q . However, here we see a possible drawback of 
assuming an exact power-law expansion, equation ( f4.12| ): these background solutions 
constrain the range of a, ocq, and uq, and from this analysis it is not obvious whether 
more general, quasi-power-law backgrounds will admit a different parameter space. 
We will leave also this question to future investigations. For the time being, we show 
that there exist quasi-power-law expansions as exact solutions and we briefly sketch 



a profile corresponding to a perturbed de Sitter background, 



— 1/r = % c and 



a q 7^ 0. Assuming that y? q = y? q o, V q = const, equation (|4.11D yields, as before, either 
V^qO = or m 2 = cr(3 — a), while equation QL10|) becomes 



1 + a 8ttG 



4nGa V c + V q 
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The solution is 



a q (t) 



(4.28) 
(4.29) 



where E is the exponential integral function. For b% > and integer a > 0, 

a e blT + Pol[0(r CT - 1 )] 



-cr+l 
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The last term, a polynomial of degree a — 1, dominates at early times (r — > — oo) 
and a q ~ |t|~ 2 , while at late times (r — >■ 0) one has a q ~ |r| _1_cr . 

When a = 0, at early times Eifor) ~ e _b2T /(&2 r ) and a q ~ |t| -2 , while at late 
times £ , i(6 2 t) ~ — In |r| and a q ~ In |t|/|t|. 

5. Scalar perturbations 

5.1 Counterterms 

As mentioned in the introduction, the corrected perturbed equations feature coun- 
terterms proportional to Sp\ in addition to the primary correction functions a and 
which guarantee consistency of the constraint algebra at 0(<5pi) order. Consistency 
in a given scheme then uniquely relates the counterterms, which all vanish classically, 
to the primary correction functions, but also restricts the range of parameters in a 
and v. Before using the counterterms in perturbation equations, we evaluate these 
consistency conditions in relation with table |l|. In the following, we shall rewrite 
the counterterms and equations of motion of according to the #pi-expansion. To 
keep notation light, background quantities will not be denoted with bars as in |J. 
Also, contrary to this reference we shall expand all intermediate expressions to linear 
order in counterterms (for instance, (1 + /)(1 + h) = 1 + / + h + 0(5^), and so on). 
Explicitly, the counterterms are 



./ 
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fi 




(5.1) 



3 dlnp 




(5.2) 
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(5.3) 



and 



1 d In a d In v 2 d 2 In v 



(5.4) 



9i = o 



3 d In p d In p 9 d In p 2 




h = h ~ 9i 



(5.5) 
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There is also the extra consistency condition 



, d/ 3 
dlnp 



+ 3(/ 3 -/) = 



which makes some of the parameters dependent: 



(5.6) 



so that 



0i 

/a 



1 3 
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1 3«o 
2a-3 



+ 1 ) - a 



^0 



6 



1 



5pi , 
5pi 



5pi , 



(5.7) 
(5.8) 
(5.9) 



the last expression being valid only if o ^ 3. 

It is interesting to notice that, for the second but not last time, the minisuper- 
space and first lattice parametrization ( |3.13| )-( |3.15| ) show an incompatibility with 
independent results: in these parametrizations, equation ( |5.6| ) is never respected. 
However, with the new range (|3.20|) for a, the equation can easily be satisfied by the 
solutions of table |l|. Let us compare case by case with equation ( |5.6|) . The solution 
in the first line of the table (equation (|4.24|) ) is also an exact solution of (|5.6|) . This is 
the limiting case of the solution in the second line (equation ( |4.25|) ) for a <C 1, giving 
tto ~ vq; for general values a < 3 (e.g., x = with m < 3 or x = — 1/4 with m < 6) 
solutions exist in this class for positive ao and u , while a® = when a = 3. The 
solutions in the third line of table [I] are already excluded by the general constraint 
«o, vq > 0. 6 Finally, the last solution in the table, equation ( f4.26| ) combined with 
(pTB]), is non-trivial and inflationary only if o = 3, but this collapses to de Sitter, 
e = 0. To summarize, our solutions will span the range 



< o < 3, 



(5.10) 



with preference to the extremum values if «o and uq are positive. In fact, e > 1/3 
is not very small for 1 < a < 2, which should lead to unviable deviations from scale 
invar iance. 



3 Even ignoring the constraint on the sign, these solutions would be inconsistent or trivial. For 



3 < a < 6, both ( 4.2 5| ) and (5ij) require ao and vq to have opposite sign, but while \vq\ 3> |ao| for 
( 4.25 ), equation ( 57q ) asks them to be of about the same magnitude. When a — 6, ao = = v^. If 
a > 6, equation ( |5.6[ ) requires ao and vq to have the same sign, in contrast with ( f.25| ); so again 
a = = vq. 
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5.2 Scalar perturbation equations 



An inhomogeneous perturbation 5(p in the scalar field induces two gauge-invariant 
scalar modes $ and \I/ in the metric, which are proportional to each other |6J: 



$ = (1 + h)ty. 



(5.11) 



After solving the equations for the expression for $ will be readily obtained via 
equation ( 5.11| ). The scalar- field perturbation and \I/ are related by the diffeomor- 
phism constraint (equation (90) of M) 



AnG-v'fy = + (i + / + nym . 

V 

Using this equation, one can show that the perturbed equation for \l/ is 7 

+ UF^' - (s 2 A + m%) = 0. 

The friction term is 
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where 



The (squared) propagation speed of the perturbation is 



s 2 = a 2 (l-/ 3 ) = l + x 5 : 



where 



3 V6 



pi 



(5.15) 
(5.16) 

(5.17) 
(5.18) 



Finally, the effective mass term is 
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(5.19) 



7 This is obtained by combining our equations (5.11) and (5.12) with equation (82) of [pi. In 
equation (82) one should correct the typographical error a 2 A<!> — > a 2 A^ [E7I: 



a 2 A* - SH(1 + /)[*' + (1 + / + hyWSi] = 4ttG-(1 + h) [p'Sip' - <^' 2 (1 + A + h)V + vpV tV> 6(p] , 

(5.13) 

where A is the Laplacian in comoving spatial coordinates. Also, in equation (91) of || one must re- 
place a 2 (l + h)A^ with a 2 A^. Unfortunately, this typo propagated in some of the other equations, 
so our results supersede those of || when in disagreement. 
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where 



= [2(e - rj) + (1 + (t)]oq + v (| - l) . 



(5.20) 



Taking equation (84) of |J, expanding it to leading order in quantum corrections, 
and making use of equations (|3.30Q , Q3.31]) , and (|5.11|) , one obtains the perturbed 
Klein-Gordon equation for the gauge-invariant perturbation 5(p: 



where 
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(5.22a) 
(5.22b) 

(5.22c) 

(5.22d) 
(5.22e) 
(5.22f) 



Before proceeding, we notice a potentially serious problem. In order to avoid 
superluminal propagation of signals, one should impose 



s < a 



(5.23) 



where we used the fact that photons propagate with speed a greater than the classical 
one M. Then, it should be / 3 > 0. For this to happen, we can have: 



< cr < 3: equation (|5.9|) imposes «o < 0. 



a = 3: equations ( |5.5| ) and ( |5.6| ) impose, respectively, u < and ao = 0. 
3 < a < 6: equation (|5.6|) imposes a and u to have opposite sign. 



<j = 6: equations ( |5.9|) and ( |5.6|) impose, respectively, a > and u = 0; this 
case is allowed. 



a > 6: this case, too, is allowed, with both oq and u strictly positive. 
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Unfortunately, for non- negative ao and u , / 3 is negative unless a be large enough, 
and this condition is hardly compatible with inflation; see table [[} (In ||47|| ) the values 
of parameters given for subluminal evolution correspond to the case a > 6 here.) 

To check whether superluminal propagation is an artifact of linear perturbation 
theory or of the expansion in 5p\, one should go beyond linear order in both expan- 
sions. The covariant formalism of non-linear perturbation theory could be a useful 
tool for analyzing the consistency of the effective constraint algebra. A possibility 
is that holonomy corrections, which we have ignored, would play an important role 
in this issue, which we shall put aside in this paper. However, even if this were the 
case in some regimes, one can always find initial conditions so as to have dominant 
inverse-volume corrections; thus, superluminal velocities might constitute a concep- 
tual problem with implications for the stability of the theory as a whole. On the 
other hand, we note that inflationary models based on superluminally propagating 
fields have been consistently formulated |HB, |{|. The case of superluminal motion 
found here therefore does not necessarily mean a severe problem. 

We reemphasize the importance of equation ( |5.6j ) and of the counterterms it 
comes from. It rules out the minisuperspace-related parametrizations and severely 
restricts the lattice one. In this way, consistency alone already subjects the theory 
to strict tests even before evaluating the phenomenology, to which we turn now. 



5.3 Conservation of curvature perturbation 

The gauge-invariant linear comoving curvature perturbation is 



1Z 



ifi' 

* + — [!- -^pi 5<p . 



(5.24) 
(5.25) 



ip' \ 6 

In the absence of counterterms, conservation of the energy-momentum tensor implies 



that 1Z is constant at large scales [41]. One may ask if this result, which is not 



obvious in Hamiltonian formalism and for equation ( 5.24j) , holds also in semi-classical 
LQC. To check it, we invert equation ( |5.13|) with respect to Sip 1 and employ (|5.12|) . 
Differentiating TZ with respect to conformal time, we obtain 



n' = {av + f-h-h 
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(5.26) 
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Here, after using ( 3.26Q , the prefactor of the parameter e (which we are not assuming 
to be small) vanishes by virtue of equation Q5.2|) , which also implies (/' — f[)/'H = 
(2/3)d 2 lnz//dlnp 2 . Together with f - f 3 = -(2/3)dlnz//dlnp+ (l/3)dlna/dlnp + 
(2/9)d 2 In u/d Inp 2 from ( |5.5|) and (|5.4|) , all terms are zero. The resulting conserva- 
tion of power is consistent with the picture of an effective perturbed FRW geometry 
that models the dynamics of a nearly isotropic universe in the presence of correc- 
tions from loop quantum gravity. Quantum-geometry corrections from this theory, 
with the perturbation equations used here, have been shown to deform not just the 
dynamics but also the underlying spacetime structure, inferred by an analysis of 
the algebra of constraints. The gauge transformations they generate no longer cor- 
respond to pure coordinate transformations because they do not obey strictly the 
classical algebra of spacetime deformations. As a consequence, classical results about 
the conservation of power may no longer apply. As seen here, the linear curvature 
perturbation is nevertheless conserved on large scales. This observation demonstrates 
that perturbations in the presence of quantum corrections can still be seen as those 
of an effective line element: The large-scale curvature perturbation of an FRW line 
element in conformal time amounts simply to a spatially constant rescaling of the 
scale factor, which should not be subject to non-trivial dynamics. By being con- 
served also in the presence of quantum corrections, the interpretation of the effective 
geometry as a line element is still meaningful. 

The result ( |5.26|) is due to the delicate cancellations between counterterms. 
Therefore, 



(5.27) 



and the curvature perturbation is conserved at large scales. 
5.4 Mukhanov equation 

Conservation of 1Z strongly suggests that one can write a simple Mukhanov equation 
in the variable 

u = zll, (5.28) 

where z is some background function. We can anticipate the main result with a 
very efficient trick, and then confirm it via a standard but tedious calculation. The 
trick is to notice that, at super-horizon scales, the co moving curvature perturbation 
is approximately constant, so that u" « z"lZ and 




u 



— u ~ . 

z 



The objective now is to find this friction-free Mukhanov equation from the perturbed 
equations of motion. Start from equation ( 5.21 ) and choose for simplicity a spatially 
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flat slice where \I> ~ 0. In order to remove the friction term, we need to define a field 
u = a(l — (35pi)5(p, where (3 = B w /a. Then, 



u 



1-135 



pi 



5<p" + 2n(l + B 10 5 P1 )S(p' + 



Comparing with the Mukhanov variable (|5.28|) and equation (|5.25|), one finds 
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(5.29) 



The only missing term in the Mukhanov equation is the Laplacian, with coefficient 
s 2 as an inspection of equation (|5.21| ) immediately shows. Thus we obtain 



u 



s 2 A + — )u = 



(5.30) 



a result valid exactly at all scales and at the linear perturbative level. The rigor- 
ous calculation begins with the Mukhanov variable ( |5.28| ) and equation ( |5.29| ) with 
unknown (3. Differentiating u twice, using equations ( |5.14p and ( |5.21| ), and using 
equation (|5.12|) to develop the term, we obtain equation (|5.30|) plus just one extra 
term: 



s 2 A 



z 
z 



u+2a'H5pi 
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«0 
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The extra term vanishes if f3 is chosen as above. 

It may seem that equation ( |5.30| ) is not covariant since only the spatial-derivative 
term is corrected. However, despite appearance this is not the case: The quantum- 
corrected equations of motion correspond to a deformed algebra of constraints as 
found in 0, and the constraints determine what form gauge transformations take. 
In general relativity, the gauge transformations are spacetime diffeomorphisms or 
changes of coordinates whose classical form implements the usual notion of covari- 
ance. With corrected constraints obeying a deformed algebra, the gauge transforma- 
tions are not of the classical form, and they do not correspond to the usual notion 
of coordinate changes. Even though the underlying structure of a 'quantum man- 
ifold' (perhaps non-commutative) on which the modified transformations could be 
interpreted as simple coordinate changes is unknown, the (generalized) covariance of 
( |5.30 ) under these deformed transformations is guaranteed by the derivation of the 
equations of motion used here from an anomaly-free set of constraints. 

It is quite remarkable that scalar perturbations are ultimately governed by such a 
simple equation as ( |5.30 ). However, the existence of one Mukhanov variable obeying 
one equation in closed form is not unexpected, as it could have been inferred by us- 
ing the Hamilton-Jacobi method for constrained Hamiltonian systems developed in 
PH , 1ST) . In particular, the reduced phase space obtained after solving the constraints 
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and factoring out their gauge flows has one local degree of freedom, parametrized by 
the curvature perturbation and its conjugate momentum. There must be a closed 
form for the dynamics on this reduced phase space, such that Hamiltonian first-order 
equations of motion exist involving only TZ and its momentum, and they are linear 
thanks to the linear perturbation scheme used. As always, first-order Hamiltonian 
equations of motion can be expressed as one second-order equation for the config- 
uration variable, here TZ. The second-order equation in general may have terms 
involving TZ", TZ' , as well as 7Z, which on large scales is an ordinary differential equa- 
tion with gradient-free coefficients (in momentum space, they are independent of the 
wave number k defined below). Thus, one may eliminate the last term involving TZ 
by substituting ylZ for TZ for a suitable background function y, and a constant mode 
for yTZ results. As a consequence, there must be a conserved quantity yTZ whose ex- 
istence can be seen without any detailed calculations. Details are required to derive 
the form of y, and the non-trivial result found here is that y = 1. 

From now on we expand linear perturbations in momentum space, a subscript 
k indicating modes with comoving wavelength 2n/k. The Laplacian becomes A — > 
—k 2 , and the Mukhanov equation 




(5.31) 



The effective mass term is a combination of slow-roll parameters and quantum cor- 
rections. In fact, 
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(5.32) 
(5.33) 



» u = ^ + u (a - 3) + + ^ - 2u ) e + 2(u - a )rj . ( r,.3 I ) 



When the slow-roll parameters are constant classically, as in any of the solutions of 
section |], one has 

zT_ _ 4^ 2 - 1 + 4/i 2 5pi 

z 



where 



/i2 = an 



4r 2 



cr(4n — an + 1)^-^ — /i„ 



(5.35) 



(5.36) 



An exact solution of the Mukhanov equation does exist but it is too complicated and 
not very instructive. We proceed to solve this equation asymptotically. 
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5.5 Asymptotic solutions 

The moment of horizon crossing is, as usual, defined when the effective mass term 
equals the Laplacian term. Up to numerical factors, this happens when 

fc|r| = l, (5.37) 

as in standard inflation. Super-horizon modes are characterized by k\r\ <C 1, while 
modes well inside the horizon have k\r\ 3> 1. At large scales, we can ignore the k 2 



term in equation ( |5.31| ), so that 



u k fc|r ^ <1 C(k)z , (5.38) 

where C(k) is a normalization constant. To determine it, we must find the asymptotic 
behaviour of u at small scales. There, one can ignore the mass term and consider 
the equation 

< + (l + xM^V~0. (5.39) 

Since all the analysis is valid only at first order in the quantum corrections, it is 
consistent to look for short-wavelength solutions of the form 

Wfc»«0) = u c (k, t) [1 + y(k, r)5 P i] , (5.40) 

where u c is the solution of the classical Mukhanov equation and y is some function. In 
particular, the only choice compatible with the Bunch-Davies vacuum in the infinite 
past is an incoming plane wave, 

Bo = (,41) 



The normalization here is the classical one, which one might have to change for a 
vacuum matter state in a quantum geometry. In particular, the correction function 
v multiplies the kinetic term of the scalar Hamiltonian, and thus affects the value of 
vacuum fluctuations. By the ansatz ( |5.40| ), all these effects will be included once the 



equation of motion for y is solved. 

Plugging the ansatz ( |5.40|) into (|5.39| ) we obtain an inhomogeneous equation for 



the function y: 

y" - 2{aU + ik)y' + 2\kaUy + xk 2 = , (5.42) 

where, for consistency, we have dropped the mass term <r'H 2 (cr + e — l)y. At this 
point we expand y in a power series, 

+oo 

y = E y^ m > ( 5 - 43 ) 

m=0 
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and we pick a power- law background, % = n/r. Then, equation (|5.42|) is 
= 2an(iky - y x )- + [2ik(an - l)yi + 2(1 - 2an)y 2 + xk 2 ] 

T 

+00 

+ ^2[2ik(crn - m)y m - (m+ \){2an - m)y m+ i]T m ~ l . 



(5.44) 



m=2 



These terms must vanish order by order separately. If a = 0, then yo is unconstrained, 
while yx = —ifa/k + k\/2) and y m = 2{2ik) m ~ 2 y2/m\ for all m > 2. Summing the 



series, one obtains 
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We can argue that y 2 — because otherwise it in equation (|5.40|) would contain also 
an outgoing mode e +lfcr . If a 7^ 0, one obtains the following conditions: 
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(5.45) 
(5.46) 

(5.47) 



The recursive relation would determine the sum of the series ( |5.43j ), but analytic 
continuation to the case a = requires y m — for m > 2. This fixes both 1/0 and 1/1 
and the result is 
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2 an - 1 



1 + ifcr) . 



(5.48) 



The normalization of equation (|5.38|) is thus obtained by imposing the junction con- 
dition |wfc>%| = |wfc<%| at horizon crossing. Then, 
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(5.49) 



where z(k) = z(r = —1/k) and Sp\(k) = 6pi(t = — 1/k) oc k T 



5.6 Scalar spectrum, spectral index and running 

The scalar spectrum is defined as the two-point correlation function of the curvature 
perturbation 1Z over a momentum ensemble at large scales, evaluated at horizon 
crossing: 



2w 2 z 2 
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(5.50) 
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The scalar spectral index is defined as 



ra„ — 1 



dlnP s 
din A; 



For a power-law background, we have 
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(5.51) 

(5.52) 
(5.53) 



We can obtain more portable expressions by writing the spectrum on a general quasi- 
de Sitter background. Since, using ( |3.26| ), 
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where we used k = H and 
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(5.54) 



(5.55) 



Notice that if a = 0, the quantum correction is constant and the only change with 
respect to the classical case is the normalization of the spectrum. In that case, 
7s = — 5«o/2 could be of either sign. If 7 S ^ 0, there is a large-scale enhancement 
of power because 5p\ ~ a~°" ~ (l/|r|) _cr ~ k~ a at horizon crossing. The magnitude of 
the effect depends on the value of 7 S but we notice that, if a ^ 0, 7 S ~ 3a /(2e) unless 
vq ~ «o/e S> «o- Therefore, one could obtain a sizable enhancement unless «o < 
O(e) (which is the case, typically). Since this enhancement is of potential interest 
for comparisons with observations, we trace back where the inverse of the slow-roll 
parameter in the expression for 7 S came from. It arises due to the e-independent 
term in z 2 above, which in turn is a direct consequence of the presence of gravity 
corrections in the Raychaudhuri equation ( |3.3| ), as opposed to just stress-energy 
modifications. As with several other key phenomena pointed out here, this feature 
is a consequence of corrections to the structure of spacetime geometry: corrections 
in the terms ~H', ~H 2 of the Raychaudhuri equation (or the isotropic Einstein tensor) 
can be obtained only by changing the geometrical form of gravity. 

Momentum derivatives are converted into conformal time derivatives via 

d 1 d 
din A; ~ Hd^' 
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so that the scalar index is 



1 = 2r] - 4e + a lna 5 
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where 
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(5.56) 



(5.57) 



Since the quantum correction is small, the scalar index does not deviate too much 
from scale invariance. If a = 0, there are no corrections at all. If a ^ 0, the sign of 
the correction depends on the choice of the parameters in the parameter space. We 
have seen that power- law/quasi de Sitter solutions have cr < 3, so it is immediate 
to associate scale invariance with small values of a. The naturalness of this range is 
further stressed in the concluding section by an independent argument. 

Interestingly, the running of the spectral index is dominated by the quantum 
correction (unless a = 0): 
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(5.58) 
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This result signals a qualitative departure from classical inflation, since the quan- 
tum correction may be larger than 0(e 2 ). The details will depend on the chosen 
background, as the slow-roll parameter themselves can contain quantum corrections. 



6. Tensor perturbations 

The linearized equation of motion for tensor modes has been computed in H and 



solved in |TT| for quasi-classical inverse volume corrections. In this section we re- 
view and improve these results, eventually obtaining the cosmological consistency 
relations. 

6.1 Mukhanov equation 

When only inverse-volume corrections are taken into account and in the absence of 
anisotropic stress, the equation of motion for the individual tensor mode h k is P 

K + 2u(l- ^) h' k + a 2 k 2 h k = . (6.1) 

Defining 

w k = ah k , a = a (l - y^Pi) , (6.2) 
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we get the Mukhanov equation 




(6.3) 



where 
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Equation Q3.3D is formally identical to the scalar Mukhanov equation and the analysis 
is exactly the same up to the substitutions 

z — y a , x ~* 2cto • 
The final result is the analogue of equation ( |5.49| ), 
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6.2 Tensor spectrum, spectral index and running 

The tensor spectrum is 

32Gk 3 



so that in de Sitter (n 
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where 
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As for the scalar spectrum, barring special values of the parameters (7t = when 
a = 1 or a = 0) there is a power enhancement at large scales because of 5p\ ~ 
k~ a , albeit the prefactor might not be as large as in equation ( |5.55|) . This type of 
enhancement has been seen in earlier numerical studies of the LQC tensor power 
spectrum, but it is difficult to exploit it observationally due to limitations by cosmic 
variance. 

The tensor index and its running are 
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(6.10) 



- 32 - 



so that 



n t = -2e - a^ t &_ 
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(6.11) 



and 



a t = -4e(e - -q) + a(2e + aj t )5 P i . 



(6.12) 



6.3 Tensor-to-scalar ratio 

The last piece of information we want to extract is the tensor-to-scalar ratio 

V, 



r = 



From equations (|5.54 ) and fl6.8|) one obtains 



r = 16e[l + ( 7 t - 7s) <*pi] 
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which yields the consistency relation 



r = -8{n t + [n t (7t - 7s) + c7t]<W 



(6.15) 



Here we implicitly assumed that 7 S is not too large, so that the expansion in <5pi is 
still meaningful. In quasi de Sitter regime e < 1, so that 7 S = 0(1) if aa ~ 0(e). 
This means that either a or ao or both should be small. 

Unless a = or o = 1 (for which 7 t = 0), the tensor-to-scalar ratio is no longer 
proportional to the tensor index. Detection of a non-zero r would require either 
a consistent deviation from de Sitter in standard cosmology or a sufficiently large 
quantum correction in de Sitter LQC. 

As already explained for ( |5.30|) , equations (|5.31|) and ( |6.3|) are covariant under the 
deformed transformations generated by the anomaly-free set of corrected constraints. 
The deformation gives rise to a new type of quantum effects which could not be 
present for higher-curvature effective actions usually expected of quantum gravity; it 
is possible only thanks to quantum corrections to the geometry of space or even the 
manifold structure. The Mukhanov equations for scalar and tensor modes are not 
only corrected, they also acquire corrections of different forms. The scalar equations 



has a correction given by s 1 



a 



^1 ~ f?i)i while the tensor equation is corrected by 



a 2 . The counterterm f% cannot typically be set to zero, and so the corrections for 
scalar and tensor modes differ. This difference, in turn, makes possible changes to 
the tensor-to-scalar ratio which may provide a key signature of loop quantum gravity. 
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7. Discussion 



To summarize the main results of this paper: 

• Different parametrizations of LQC inverse- volume quantum corrections are pre- 
sented. A particular case is the so-called mini-superspace parametrization. 
Lattice refinement parametrizations are introduced in a self-contained way re- 
organizing a number of related results in the literature. 

• Cosmological solutions with power-law scale factor are found for the LQC back- 
ground equations of motion in the presence of small inverse- volume corrections 
Sp\. Their form depends on the parameter range, they generalize the unique 
classical power-law solution with exponential potential and are first order in 
the quantum corrections. These solutions are further restricted by an anomaly 
cancellation condition and the requirement of slow roll. 

• For the first time, the already-known equations of motion for scalar perturba- 
tions with inverse-volume corrections are consistently rewritten as first-order 
expressions in dpi and combined into a single Mukhanov equation. The LQC 
comoving curvature perturbation is shown to be conserved at large scales, just 
as its classical counterpart. 

• The Mukhanov equation is solved asymptotically and the scalar spectrum and 
index are constructed. 

• The already-known LQC Mukhanov equation for tensor perturbations is con- 
veniently rewritten and solved using the same techniques. Tensor observables 
are extracted and combined with scalar observables, thus providing the full set 
of inflationary observables in linear perturbation theory. 

As long as the slow-roll approximation is valid, the structure of the cosmological 
observables is valid for any background, although the coefficients of the quantum 
corrections themselves do depend on the background. In this final section we discuss 
how they can be used to restrict models of loop quantum cosmology, making the 



framework falsifiable. Details will be provided in separate publications p2\ ■ For such 
an endeavor, it is crucial to obtain independent information on the main correction 
parameter Sp\ and on different versions of the parametrization. For instance, as seen 



in section p. 1.3 , a combination with holonomy corrections is interesting and shows a 



powerful interplay between these main two types of quantum-geometry corrections. 

If inflation is assumed, the minisuperspace parametrization is under tight pres- 
sure: neither consistent power-law background solutions nor a nearly scale-free spec- 
trum can be found in that case (unless 5p\ be very small; see below). On the other 
hand, the minisuperspace parametrization may still be viable if an alternative sce- 
nario of structure formation can be found. In this context, holonomy corrections are 
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of particular interest not only by providing an additional consistency condition in 
combination with inverse-volume corrections, but also because they can easily trig- 
ger bounces at least in isotropic models whose matter energy is dominated by the 
kinetic term. (In general, it has not been shown that isotropic bounces occur as a 
natural consequence of holonomy corrections.) It would thus be of interest to de- 
velop linear perturbation equations around those models and analyze the structure 
evolution through the bounce, or perhaps new scenarios providing the generation 
of structure during a phase before (not after) the big bang. However, compared to 
inverse-triad corrections such ideas are currently hampered by several major diffi- 
culties: (i) Holonomy corrections have so far not been implemented in consistent 
deformations of linear perturbation equations, (ii) Strong quantum-geometry cor- 
rections are required to evolve through the bounce; no expansion in parameters such 
as Sp\ used here could be done, (iii) There are several indications as to the strong 



sensitivity of evolution through the bounce to initial conditions of perturbations [[53 
or even the quantum state J51] , discussed in the context of cosmic forgetfulness. 
The lattice parametrization is consistent with inflation and holonomy corrections 



at the homogeneous level and could yield strong effects according to section p. 1.3 



because a can be small, but that again depends on the details of Spy. We have seen 
that <5pi is an eigenfunction of the operator d/d In k with eigenvalue —a, so observables 
of higher order in the slow- roll parameters (e.g., the index running) are corrected by 
a term which is always of the form 0(a n )Sp\: it is first-order in Sp\ and n-th order in 
a. If a = 0(1), this quantum correction is equally important at any slow-roll order, 
if not increasing with the order. This situation does not seem natural inasmuch as it 
would imply that higher-order k derivatives of the inflationary spectra are all on the 
same footing. Then, the notion that the spectra can be approximated by a power law 
would have to be abandoned. On the other hand, if o <C 1 the quantum correction 
is suppressed by higher and higher powers of cr, so there is a sort of balancing effect 
which keeps 0(a n )Sp\ small at all orders in the slow- roll parameters. This leads to 
the speculation that small values of a are more sensible, because for large a quantum 
corrections would dominate in higher-order observables. 

With cr <c 1 preferred, inverse-volume corrections are of the form 1 + cp~ a / 2 
with a small exponent a. They affect not only the equations for an expanding uni- 
verse but also the dispersion relations of waves propagating in a quantum spacetime. 
Corrections for these equations are more difficult to derive because the situation is 
not as symmetric as the one of perturbations of an isotropic spacetime. But if they 
turned out to be of a similar form (Sp\ with small exponent a) for a variable related 
to the particle's energy, severe observational pressure could be put on loop quantum 
gravity by a combination of cosmological and astroparticle observations. 

The first-order cosmological observables already give a wealth of information 
about the early universe but it is natural to ask oneself what happens at second order, 
e.g., when looking at the bispectrum and possible non-Gaussian signatures. Posed in 
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LQC, this question is less harmless than in classical general relativity. In fact, when 
going to higher orders in perturbation theory one would get more parameters initially 
because there are more options for counterterms. The counterterms would then be 
fixed by consistency conditions. This would not necessarily add extra conditions for 
the parameters arising at lower orders, but something like this might happen. We 
do not know yet whether loop quantum gravity as a whole is consistent, so at some 
point parameters might be overconstrained. If that were the case, at least in LQC, 
one would have to understand if a non-perturbative cancellation of anomalies (which 
one did not see at the perturbative level) takes place. Therefore, an extension of our 
results to a second-order analysis would be most welcome not only for the purpose 
of finding the bispectrum, but also in order to further check the self-consistency of 
the theory. 

Before we conclude, we emphasize that we have considered in detail only one type 
of corrections (inverse- volume) and no complete set of effective equations implement- 
ing all the effects expected from loop quantum gravity. Even so, the conclusions we 
draw are reliable because they point out characteristic phenomena from the correc- 
tions considered. An elimination of these effects by including other phenomena (most 
importantly, those due to holonomies) can be expected only under very fine-tuned 
conditions. The equations provided here can thus be used to place bounds on the 
free parameters of loop quantum gravity, and to rule out some parametrizations as 
extensively done in this paper. In all these cases, we see the importance of sufficiently 
general parametrizations for models to be able to stand up to phenomenological pres- 
sure; conceptual preferences or 'natural' choices of parameters may not always be 
the ones that survive stringent analysis. As our examples demonstrate, it is impor- 
tant to combine constraints from different sources. For instance, the minisuperspace 
parametrization is ruled out by (i) the consistency condition provided by anomaly- 
freedom, (ii) the interplay of inverse- volume with holonomy corrections, and (iii) the 
phenomenological requirement of a nearly scale-invariant spectrum. A single incon- 
sistency could always be evaded by questioning the condition violated, especially in 
a situation in which no tight derivations from an underlying full theory exist. But as 
inconsistencies pile up, models should eventually be dropped. In this way, models of 
loop quantum gravity and loop quantum cosmology are already falsifiable not just 
by internal consistency considerations but also by comparison with observations. 

Our results are interesting also because they highlighted a number of issues which 
are definitely worthy of further attention: 

• The consistency of the slow-roll background solutions, of the anomaly cancel- 
lations, and of the physical observables only with respect to the second lattice 
parametrization urges us to study the latter in greater detail. Such an en- 
deavor requires a better understanding of the full theory and its reduction to 
perturbations around isotropic models. 
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• Perturbations can propagate with superluminal speed. Either this is an artifact 
of linear perturbation theory, is curable with a particular choice of the param- 
eters, or may give rise to non-standard inflationary scenarios as in ||49|| . The 
minisuperspace parametrization is safe as / 3 > in that case (large a > 3), 
while the lattice refinement parametrization with o <C 1 requires «o < 0. This 
observation may be the one putting the most severe constraints on the lattice 
refinement parametrization, while the minisuperspace parametrization is under 
much stronger pressure from other consistency conditions. 

We believe that addressing these points and the mutual tension between different 
parametrizations and the physical viability of the perturbations will stimulate the 
advance in the field. 
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